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Dual Opposed Stirling Convertors 

• High Efficiency, Low Mass Space Power 





Simulate Engine Only 
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Operating Conditions 
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Hot-End Temperature, K 

923 

Cold-End Temperature, K 

353 

Ambient Temperature, K 

293 

Frequency, Hz 

80 

Mean Pressure, Pa 

2.429E+6 

Power Piston Amplitude, mm 








Whole-Engine Simulation 
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PV Power, W 


Heat In, W 


Heat Out, W 


PV Efficiency 


Pressure Ratio 


Axisymmetric 

Simulation 
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Regen. AP., Pa 


Heater P, Pa 


Cooler P, Pa 


Pressure Amp., Pa 225269 


Mean Pressure, Pa 2429130 
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Comparison with Experiment 
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Simulation 
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Pressure and Heat Transfer over Cycle 




Expansion and Compression Pressure 


9 




Regenerator Pressure Drop 
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Solved Whole-Engine On 32 Processor Cluster 
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Multidimensional Analysis Tools 

• Modified CAST, CFD-ACE, Fluent, STAR- 
HPC, CFX/ANSYS, others ... 

• All low order SIMPLE/PISO based 


• New high-order codes being developed 




Third Order Accuracy 
Fifth Order Accuracy 
Seventh Order Accuracy ./ 
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Multidimensionally Enabled Capabilities 
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• Structural, thermal, fluid, electromagnetics, 
and coupled physics. 

• 1 hour per cycle axisymmetric simulation 

• Seal & appendix gap phenomena 

• HX end effects 

• Effect of vortices in expansion & compression 

• Flexure temperatures for reliability analysis 

• Effects of slight geometrical variations 
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Cool End Heatino 
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Radial Heat Transfer Effects 
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Effects of Vortices in Expansion 
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Simulation Approach 

• Whole-engine for physical consistency 

• REV Regenerator Modeling 

• Grid Layering, Smoothness, and Quality 

• Conjugate heat transfer method adjustment 

• High-speed Low Cost Parallel Cluster 

• Lot’s of debugging.... 
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Dangers of Component Modeling 
Axially Oscillating Boundaries 
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Radially Oscillating Boundaries 
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Heat Exchanaer Oscillatina Boundar 
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Manifold & Aooendix Gao O.B 
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Cooler Manifold O.B 
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Seal Region O.B 
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Compression Manifold O.B 
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Regenerator Modeling 
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• Geometrical Shapes Affect Complexity 

• High Efficiency Regenerators Permit 
Thermal-Equilibrium Assumption/Usage 

• Add source terms to governing equations 

• Darcy-Forcheimer equation is used on 
momentum, Coefficients from experiment 

• Single Energy equation (fluid/solid averaged) 
with averaged conductivity and energy 
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Regenerator Modeling in More Detail Later 
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Gridding Approach 

• Scale geometry to compensate for 32-bit 

• Avoid tetrahedrals, use unstructured quads 

• Match cell sizes across sliding regions 

• Smoothly expand/shrink cell sizes throughout 

• Use layering, avoid adaption/remeshing 

• Utilize moving and double-sided boundary layers 

• Include as much geometry as practical 

• Clean up geometry in CAD, avoid virtual geometry 

• Over-converge the steady solution to improve grid 
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Whole Engine (700K cells) 
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Left Half Grid 
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Right Half Grid 
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Both Seal Reaions 
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Appendix Gap Grid 
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Piston Region 





Conjugate Heat Transfer Strategy 
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• Utilize coupled solver in double precision 

• Maximize energy underrelaxation factor 

• Double-sided boundary layers 

• Turn off secondary temperature gradients 

• Emphasize the coarse multigrid to avoid roundoff error 

• Utilize explicit underrelaxation of temperature 

• Maintain small time-steps (160/cycle) 

• Go to higher space and time accuracy when possible 





Rear View 


Higher Solution Speeds Achievable 

•Infiniband/High Capacity Switch 

•260 processors in a single 7 foot 
tower 

•Higher-Order Techniques 
•Quad/Octa Opteron Processors 


■256 bit computing 
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Conclusions 
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• First U.S. fully converged axisymmetric 
simulation of actual Stirling engine 

• Comparison with experiment shows less than 
1 % error on power and efficiency 

• More testing required to confirm 

• Extending to three-dimensions for full part 
testing and integration 
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Disolacer Seal Close-u 
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